Quantum measurement bounds beyond the uncertainty relations 
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We give a bound to the precision in tlie estimation of a parameter in terms of the expectation 
value of an observable. It is an extension of the Cramer-Rao inequahty and of the Heisenberg 
uncertainty relation, where the estimation precision is typicaUy bounded in terms of the variance of 
an observable. 



Quantum measurements are limited by bounds such as 
the Heisenberg uncertainty relations [J, 13] or the quan- 
tum Cramer-Rao inequality which typically con- 
strain the ability in recovering a target quantity (e.g. a 
relative phase) through the standard deviation of a con- 
jugate one (e.g. the energy) evaluated on the state of the 
probing system. Here we give a new bound related to 
the expectation value: wc show that the precision in the 
quantity cannot scale better than the inverse of the ex- 
pectation value (above a "ground state") of its conjugate 
counterpart. It is especially relevant in the expanding 
field of quantum metrology 0: it settles in the posi- 
tive the longstanding conjecture of quantum optics 0- 
, recently challenged , that the ultimate phase- 
precision limit in interfcrometry is lower bounded by the 
inverse of the total number of photons employed in the 
estimation process. 

The aim of Quantum Parameter Estimation is 
to recover the unknown value a; of a parameter that is 
written into the state of a probe system through some 
known encoding mechanism Ux. For example, we can re- 
cover the relative optical delay x among the two arms of 
a Mach-Zehnder interferometer described by its unitary 
evolution Ux using as probe a light beam fed into the in- 
terferometer. The statistical nature of quantum mechan- 
ics induces fluctuations that limit the ultimate precision 
which can be achieved (although wc can exploit quantum 
"tricks" such as entanglement and squeezing in optimiz- 
ing the st ate p reparation of the probe and/or the detec- 
tion stage [l7|). In particular, if the encoding stage is re- 
peated several times using ly identical copies of the same 
probe input state px, the root mean square error (RMSE) 
AX of the resulting estimationprocess is limited by the 
quantum Cramer- Rao bound [SHal AX ^ \ j ^ vQ{x) , 
where Q{x) is the quantum Fisher information. For 
pure probe states and unitary encoding mechanism l]x^ 
Q(x) is equal to the variance {AHY (calculated on the 
probe state) of the generator H of the transformation 
Ux = e~*^^. In this case, the Cramer- Rao bound takes 
the form 



AX ^ l/{y^AH) 



(1) 




FIG. 1: Lower bounds to the precision estimation AX as a 
function of the experimental repetitions i/. The green area in 
the graph represents the forbidden values due to our bound 
((2]). The blue (dashed- hne) area represents the forbidden val- 
ues due to the Cramer-Rao bound, or the Heisenberg uncer- 
tainty. Possible estimation strategies have precision AX 
that cannot penetrate in the colored regions. For large u the 
Cramer- Rao bound (which scales as 1/ y^) is stronger, as ex- 
pected since in this regime it is achievable. Our bound is not 
achievable in general, so that the green area may be expanded 
when considering specihc estimation strategies. [Here we used 
(H) -Eo = 0.1 (a.u.) and AH ^ 4 (a.u.).] 



to the Heisenberg uncertainty relation for conjugate vari- 
ables [H . This bound is asymptotically achievable in 
the limit of i/ ^ oo HH. 

Here we will derive a bound in terms of the expectation 
value of H, which (in the simple case of constant AX) 
takes the form (see Fig. [T]) 



AX ^ K/[iy{{H) - Eo) 



(2) 



of an uncertainty relation 0,0]. In fact, if the parameter 
X can be connected to an observable, Eq. (jlj corresponds 



where Eq is the value of a "ground state" , the minimum 
eigenvalue of H whose eigenvector is populated in the 
probe state (e.g. the ground state energy when H is the 
probe's Hamiltonian) , and k ~ 0.091 is a constant of 
order one. Our bound holds both for biased and unbiased 
measurement procedures, and for pure and mixed probe 
states. When AX is dependent on x, a constraint of the 
form (2)) can be placed on the average value of AX{x) 
evaluated on any two values x and x' of the parameter 



2 



which are sufSciently separated, namely 



AX{x) + AX{x') 



(3) 



Hence, we cannot exclude that strategies whose error AX 
depend on x may have a "sweet spot" where the bound 
(21) may be beaten Q , but inequality ^ shows that the 
average value of AX is subject to the bound. Thus, these 
strategies are of no practical use, since the sweet spot 
depends on the unknown parameter x to be estimated 
and the extremely good precision in the sweet spot must 
be counterbalanced by a correspondingly bad precision 
nearby. 

Proving the bound ([2]) in full generality is clearly not 
a trivial task since no definite relation can be established 
between v{{H) — Eq) and the term ^/vAH on which the 
Cramer-Rao bound is based. In particular, scaling ar- 
guments on v cannot be used since, on one hand, the 
value of V for which Eq. ([1]) saturates is not known (ex- 
cept in the case in which the estimation strategy is fixed 
Q, which has little fundamental relevance) and, on the 
other hand; input probe states p whose expectation val- 
ues {H) depend explicitly on v may be employed, e.g. see 
Ref. 1^. To circumvent these problems our proof is 
based on the quantum speed limit [l^l, a generaliza- 
tion of the Margolus-Levitin [l^ and Bhattacharyya 



bounds [20|, |2lJ which links the fidelity F between the 



two joint states p®" and p'^1^ to the difference x' — x oi the 
parameters x and x' imprinted on the states through the 
mapping Ux = e~'^^ [The fidelity between two states p 
and (T is defined as F = {Tr[y^^7pCT^7p]}^. A connection 
between quantum metrology and the Margolus-Levitin 
theorem was proposed in j22| . but this claim was subse- 
quently retracted in [l^.l In the case of interest here, 

the quantum speed limit |18j implies 



x\ ^ — max 
' 2 



a{F) 



u{{H) - Eo) ' ^AH 



(4) 



where the f and factors at the denominators arise 
from the fact that here we are considering v copies of 
the probe states px and px', and where a{F) ~ (i'^[F) = 
4 arccos (/F)/7r2 are the functions plotted in Fig. [2] of 
the supplementary material. The inequality (|4]) tells us 
that the parameter difference jx' — a;| induced by a trans- 
formation e~*^^ -x)H employs resources {H) — Eq 
and AH cannot be arbitrarily small (when the parameter 
X coincides with the evolution time, this sets a limit to 
the "speed" of the evolution, the quantum speed limit). 

We now give the main ideas of the proof of ([2]) by focus- 
ing on a simplified scenario, assuming pure probe states 
IV's) = Ux\tp), and unbiased estimation strategies con- 
structed in terms of projective measurements with RSME 
AX that do not depend on x (all these assumptions are 
dropped in the supplementary material). For unbiased 
estimation, x = J2i ^^'^ the RMSE coincides 



with the variance of the distribution Pj{x), i.e. AX = 
Pjix)[xj - where Pj{x) = \{xj\Tpx)'^''\^ is the 
probability of obtaining the result Xj while measuring the 
joint state \tpx)'^'^ with a projective measurement on the 
joint basis \xj). Let us consider two values x and x' of 
the parameter that are further apart than the measure- 
ment's RMSE, i.e. x'-x = 2\AX with A > 1. If no such 
X and x' exist, the estimation is extremely poor: basi- 
cally the whole domain of the parameter is smaller than 
the RMSE. Hence, for estimation strategies that are suf- 
ficiently accurate to be of interest, we can always assume 
that such a choice is possible (see below). The Tcheby- 
chev inequality states that for an arbitrary probability 
distribution p, the probability that a result x lies more 
than XAX away from the average n is upper bounded 
by namely p{\x — fi] ^ XAX) < It implies 

that the probability that measuring \'^x') ■= IV'x')®'^ the 
outcome Xj lies within XAX of the mean value associ- 
ated with x) ■= \ipx)'^'^ cannot be larger 1/A^. By the 
same reasoning, the probability that measuring l^'^;) the 
outcome Xj will lie within XAX of the mean value associ- 
ated with I'^x') cannot be larger 1/A^. This implies that 
the overlap between the states \'^x) and \'9x') cannot be 
too large: more precisely, F = |(*a:|vl'2./)|2 ^ 4/A^. Re- 
placing this expression into (j4|) (exploiting the fact that 
a and /3 are decreasing functions) we obtain 



2 A AX ^ - max 



a(4/A^) /3(4/A^) 
u{{H) - Eo) ' V^AH 



(5) 



whence we obtain ^ by optimizing over A the first term 
of the max, i.e. choosing k = sup;^ tt a(4/A^)/(4A) ~ 
0.091. The second term of the max gives rise to a 
quantum Cramer-Rao type uncertainty relation (or a 
Heisenberg uncertainty relation) which, consistently with 
the optimality of Eq. ^ ioi i/ ^ 1. has a pre- factor 
7r;3(4/A^)/(4A) which is smaller than 1 for all A. This 
means that for large ly the bound ([2]) will be asymptot- 
ically superseded by the Cramer- Rao part, which scales 
as (X l/v^ and is achievable in this regime. 

Analogous results can be obtained (see supplementary 
material) when considering more general scenarios where 
the input states of the probes are not pure, the estima- 
tion process is biased, and it is performed with arbitrary 
POVM measurements. (In the case of biased measure- 
ments, the constant k in ([2]) and ([3]) must be replaced by 
K = sup^ 7ra(4/A2)/[4(A + 1)] ~ 0.074, where a +1 term 
appears in the denominator.) In this generalized context, 
whenever the RMSE depends explicitly on the value x of 
the parameter, the result © derived above is replaced 
by the weaker relation ([3]). Such inequality clearly does 
not necessarily exclude the possibility that at a "sweet 
spot" the estimation might violate the scaling ([2]) . How- 
ever, Eq. ([3]) is still sufficient strong to exclude accura- 
cies of the form AX{x) = 1/ R(x^v{H)) where, as in 
Refs. [13, Hi, RiXyz) is a function of z which, for all x, 
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increases more than linearly, i.e. lim2_j.oo z/R{x, z) = 0. 

The bound ^ has been derived under the explicit as- 
sumption that X and x' exists such that x' — x ^ 2AAX 
for some A > 1, which requires one to have x' — x ^ 2 AX. 
This means that the estimation strategy must be good 
enough: the probe is sufficiently sensitive to the trans- 
formation Ux that it is shifted by more than AX during 
the interaction. The existence of pathological estima- 
tion strategics which violate such condition cannot be 
excluded a priori. Indeed trivial examples of this sort can 
be easily constructed, a fact which may explain the com- 
plicated history of the Heisen berg bound with claims 0- 
[l3l | and counterclaims 14-3, Ij]- 1^ should be stressed 
however, that the assumption x' — x ^ 2 AX is always 
satisfied except for extremely poor estimation strategies 
with such large errors as to be practically useless. One 
may think of repeating such a poor estimation strategy 
V > 1 times and of performing a statistical average to 
decrease its error. However, for sufficiently large v the 
error will decrease to the point in which the v repetitions 
of the poor strategy are, collectively, a good strategy, and 
hence again subject to our bounds ([2]) and ([3]). 

Our findings are particularly relevant in the field 
of quantum optics, where a controversial and longly 
debated problem [S-l^ 24 1 is to determine the scaling 
of the ultimate limit in the interferometric precision 
of estimating a phase as a function of the energy {H) 
devoted to preparing the v copies of the probes: it 
has been conjectured that the phase RMSE is 

lower bounded by the inverse of the total number of 
photons employed in the experiment, the "Heisenberg 
bound" for interferometry^. Its achievability has been 
recently proved [25j, and, in the context of quantum 
parameter estimation, it corresponds to an equation of 
the form of Eq. ([5]), choosing x = (p (the relative phase 
between the modes in the interferometer) and H = a^a 
(the number operator). The validity of this bound has 
been questioned several times Il4l-ll6|. l24| . In particular 
schemes have been proposed [ij, |2J] that apparently 
permit better scalings in the achievable RMSE (for 
instance AX sa {v{H))~'^ with 7 > 1). None of these 
protocols have conclusively proved such scalings for 
arbitrary values of the parameter x, but a sound, clear 
argument against the possibility of breaking the 7 = 1 
scaling of Eq. ([2]) was missing up to now. Our results 
validate the Heisenberg bound by showing that it applies 
to all those estimation strategies whose RMSE AX do 
not depend on the value of the parameter x, and that 
the remaining strategies can only have good precision 



^ This "Heisenberg" bound [a 4l3ll should not be confused with 
the Heisenberg scaling defined for general quantum estimation 
problem 0| in which the y/u at the denominator of Eq. |[TJ is 
replaced by v by feeding the v inputs with entangled input states 



for isolated (hence practically useless) values of the 
unknown parameter x. 
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Supplementary material 

Our bound refers to the estimation of the value x of 
a real parameter X that identifies a unitary transforma- 
tion Ux ~ e~'^^, generated by an Hermitian operator H . 
The usual setting in quantum channel parameter estima- 
tion (see for a recent review) is to prepare v copies of 
a probe system in a fiducial state p, apply the mapping 
Ux to each of them as p ^ px = UxpUx\ and then per- 
form a (possibly joint) measurement on the joint output 
state p®'^ , the measurement being described by a generic 
Positive Operator- Valued Measure (POVM) of elements 
{Ej}. [The possibility of applying a joint transformation 
on the u probes before the interaction Ux (e.g. to entan- 
gle them as studied in [17[) can also be considered, but it 
is useless in this context, since it will not increase the lin- 
ear scaling in u of the term v{{H) — Eq) that governs our 
bounds.] The result j of the measurement is finally used 
to recover the quantity x through some data processing 
which assigns to each outcome j of the POVM a value 
Xj which represents the estimation of x. The accuracy of 
the process can be gauged by the RMSE of the problem, 
i.e. by the quantity 



AX := ^ V P,{x)\x, - x]2 = ^S^X + {x-xy, (6) 



where Pj{x) — Ti[Ejpf'''] is the probability of getting the 
outcome j when measuring pf"', x :~ Pj{x)xj is the 
average of the estimator function, and where 



S^X :=^P,-(x)[x,-x]2 



(7) 



e.g. see Ref. 
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is the variance of the random variable Xj . The estima- 
tion is said to be unbiased if x coincides with the real 
value X, i.e. x = x, so that, in this case, AX coincides 
with SX. General estimators however may be biased with 
a; 7^ a;, so that AX > 6X (in this case, they are called 
asymptotically unbiased if x converges to x in the limit 
v — > 00). 

In the main text we restricted our analysis to pure 
states of the probe p = and focused on projective 

measurements associated to unbiased estimation proce- 
dures whose RMSE AX is independent on x. Here we 
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extend the proof to drop the above simphfying assump- 
tions, considering a generic (non necessarily unbiased) es- 
timation process which allows one to determine the value 
of the real parameter X associated with the non neces- 
sarily pure input state p. 

Take two values x and x' of X such that their associ- 
ated RMSE verifies the following constraints 

AX{x) =^ , (8) 
\x-x'\ = {X + l)[AX{x)+AX{x')], (9) 

for some fixed value A greater than 1 (the right hand 
side of Eq. ^ can be replaced by \[AX{x) + AX{x')] if 
the estimation is unbiased). In these expressions lS.X{x) 
and AX (a;') are the RMSE of the estimation evaluated 
through Eq. ([5]) on the output states p^" and p®!" re- 
spectively (to include the most general scenario we do 
allow them to depend explicitly on the values taken by 
the parameter X). In the case in which the estimation is 
asymptotically unbiased and the quantum Fisher infor- 
mation Q{x) of the problem takes finite values, the con- 
dition ^ is always guaranteed by the quantum Cramer- 



Rao bound [3|-|6| (but notice that our proof holds also 
if the quantum Cramer-Rao bound does not apply - in 
particular, we do not require the estimation to be asymp- 
totically unbiased). The condition ([9]) on the other hand 
is verified by any estimation procedure which achieves a 
reasonable level of accuracy: indeed, if it is not verified, 
then this implies that the interval over which X can span 
is not larger than twice the average RMSE achievable in 
the estimation. 

Since the fidelity between two quantum states is the 
minimum of the classical fidelity of the probability dis- 
tributions from arbitrary POVMs [2^, we can bound the 
fidelity between pf" and pf," as follows 



F := 




E 

3 



(10) 



with Pj{x) = TilEjp^"] and Pj{x') = TT[Ejpf/]. The 
right-hand-side of this expression can be bound as 



^ JT.p^(^) E Pri^') + JT.P^(^) E Pri^') 



^ jT.P^(^) + jT.Pri^')^ (11) 

where / is a subset of the domain of possible outcomes j that we will specify later, and where we used the Cauchy- 
Schwarz inequality and the fact that Pj' (^') ^ 1 ^^^^ Sj^/ Pji^) ^ 1 independently from /. Now, take / to be 

the domain of the outcomes j such that 

\xj-x'\^\6X', (12) 

where A is a positive parameter (here x' and {6X'Y are the average and the variance value of Xj computed with the 
probability distribution Pj{x')). From the Tchebychev inequality it then follows that 

^Fy(a:')<l/A^ (13) 

which gives a significant bound only when A > 1. To bound the other term on the rhs of Eq. (|lll) we notice that 
|a; — x'l ^ |x — -I- |a;' — .t'I + |a; — .t'I and use Eq. © and ^ to write 

\x-x'\ ^ {\ + l){AX + AX') -\x~ x\~\x' ~ x'\ 



= (A + l)iAX + AX') - \/A^X - 6^X - x/A^X'^S^X' > A(AX + AX') . (14) 
From Eq. (|12p we also notice that for j £ I we have 

\x — x'\ ^ \x — Xj\ + \xj — x'\ ^ |a; — Xj| -I- X6X' , (15) 
which with the previous expression gives us 

\x -xA ^ A(AX + AX') - XSX' ^ XAX ^ XSX , (16) 
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and hence (using again the Tchebychev inequality) 



(17) 



Replacing and (HT)) into (dU]) and (HII) we obtain 

F 4/A2 . (18) 
We can now employ the quantum speed limit inequality (|4]) from [isj and merge it with the condition (|9]) to obtain 



(A + 1)(AX + AX')^|.'-.| ^ - max | . ^ ^ max | ^^^^^^^ , ^ | 



, (19) 




FIG. 2: Plot of the functions Q(i*") and /3(F) appearing in 
Eq. (gl). 



where, as in the main text, we used the fact that a and 13 
are decreasing functions of their arguments, and the fact 
that the expectation and variances of H over the family 
Px is independent of x (since H is independent of x) . The 
first term of Eq. together with the first part of the 
max implies Eq. ([3|), choosing k = sup;^ 7ra(4/A^)/[4(A + 
1)] ~ 0.074, which for unbiased estimation can be re- 
placed by K = sup^ 7ra(4/A^)/[4A] ~ 0.091. In the case 
in which AX{x) = AX{x') = AX we then immediately 
obtain the bound 
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